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Using the BRM theory developed recently by Fyodorov and Mirlin we calculate the density- 
density correlators for Banded Random Matrix of infinite size. Within the accuracy of 1/6^ {b is the 
matrix bandwidth) it appears to be the same in both cases of the orthogonal and unitary symmetry. 
Moreover, its form coincides exactly with the formula obtained long ago by Gogolin for the electron 
density-density correlator in strictly ID disordered metals. 

In addition to the "fixed energy" density-density correlator considered in the solid state physics 
we calculate also the "time averaged" one, which has different properties at small separations. Our 
, predictions are compared with the existing numerical data. 

a^ 

I. INTRODUCTION 

O 

^ • The Band Random Matrix(BRM) model is a good basisQ for understanding basic features of quantum chaotic 
' systems. It looks even more attractive in view of its analytical solution obtained by Fyodorov and Mirlin in the 
framework of the Efetov's supersymmetry approach |^ . Being inspired by numerous earlier attempts to understand 
scaling properties of the BRM model, the solution is proved to be an adequate tool to address the problem. 

Recently, further extensive numerical studies of the RBM model for the real symmetric matrices (we refer to this 
case as that of orthogonal ensemble) has been done in the work Among many of interesting results concerning 
_ : the scaling properties of the BRM model, very detailed data were presented for the time- averaged density-density 
' correlator, known as the "steady state distribution" 
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■ K.t(xi,X2)^ Yuii^-^ dt\{x2,t\xi,Q)\^ =^\ti)l{x2)lpa{xi)\^, (1) 
ON ' 

P5 : where ipa are eigenfunctions of the Hamiltonian considered to be a member of BRM ensemble: Hipa = Eatpa and Ea 

J>~^' is the corresponding eigenvalue. The scaled coordinates enetering Eq.(|^) are xi^2 = ni^2/b^, with ni_2 being labels of 

, the components of BRM eigenvectors and b being the band widths. The bar stands for the averaging over the random 

Q ' matrix ensemble. 

d In @ the data for the correlator (^ has been fitted by the theoretical expression obtained long time ago by Gogolin 
^ for strictly ID Anderson model. The Gogolin's expression contains the localisation length / as a parameter. A 

J-^ reasonable agreement has been found for the length parameter value Im w 0.296^. This value is rather close to the 

I> , localization length loo\fT^_n\ = k^^ expected for the orthogonal BRM ensemble with the matrix size ^ 5^ ^ 1 n. 

However, the fit done in |^ is in fact not very legitimate. The reason is that the Gogolin expression is not relevant 

; I for the correlator and represents another quantity: a "fixed energy" correlator defined as 

a 



ICE{xi,X2;E)^J2\r^{x2,E^)^4xuE^)f6iE-E^). (2) 

a 

It is easy to see, that time-averaged correlator can be obtained from the "fixed energy" one by integration over 
the energy. Since the localization length parameter depends essentialy on the energy, the integration can modify 
the correlator form. Indeed, as we show below the Gogolin correlator with energy dependent localization parameter 
integrated over the energy is a better fit for the data. 



*e-mail address zhirov@inp.nsk.su 

^See, e.g., as well as a recent very extensive review and references therein. 

■^Note, that in Ref. |^ it was erroneously concluded that that Im differs substantially from Zoo- This is however a consequence 
of incorrect use of the localization length loo defined for the unitary (see Eq.(4),(5) in Ref. [^) instead of the orthogonal BRM 
ensemble. The localisation length Zoo for the unitary class is known to be twice as large in comparison with that for the 
orthogonal class.) 
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Actualy, the best way is to fit the data by theoretical predictions following from the BRM model itself, rather than 
to use Gogolin's formula derived formally for diferent case. Unfortunately, the relevant quantity was not calculated 
up to now. In particular, in |^ the correlator /C^ (xi , a;2 ; -E) has been calculated at points n(xi) = 0, n{x2) — N {N 
is matrix size), while to compare with data one needs such a correlator for finite separations X2 ~ xi but infinite 
matrix size N oo. 

This work is intended to fill this gap. We calculate the density-density correlator for infinite banded random 
matrices of two symmetry classes: the orthogonal ensemble, studied in ||^, and the unitary one, consisting of complex 
Hermitian banded matrices. The paper is organized as follows: after a brief overview of mapping of the BRM model 
onto the supersymmetrical cr-model in Sec.||, we consider step by step the transfer matrix solution of the cr-model 



in Sec. Ill and calculate the fixed energy correlator in Sec. IV for both symmetry cases simultaneously. As a result, 
we find that within an accuracy of 0{l/b^) both correlators coincide in their functional form. At the same time, 
the localization length is as expected twice as smaller in the orthogonal ensemble as compared to the unitary one. 
Integrating in Sec.^ the "fixed energy" correlator over the energy we complete the paper with BRM predictions for 
the time averaged correlator, the "steady state distribution". 

II. BASIC DEFINITIONS OF THE BRM PROBLEM AND ITS MAPPING ONTO THE 
SUPERSYMMETRIC NONLINEAR a-MODEL. 

Let us consider the ensemble oi N x N {N ^ 1) banded random matrices H, whose elements are independently 
distributed according to the Gaussian law with zero mean value and variances 



\H,,f ^J,,^fi\i-j\), (3) 

The band profile function /(r) vanishes exponentially or faster outside the band |« — j| < 6, 6 3> 1. In this case, as it 
was known from |^ predictions are insensitive to the particular band profile form. 

Below we consider simultaneously two symmetry classes of matrices H: the orthogonal ensemble of real symmeric 
matrices, and the unitary ensemble of hermitian ones. The former corresponds to systems which are time-reversal 
invariant , while the latter corresponds to broken invariance. 

One can get a solution of the BRM problem by mapping it onto supcrsymmetric cr-model studied in many detailes 
in earlier Efetov's works [p^|6|. Let us summarize the results of such a mapping for both symmetry classes. The 
corresponding cr-model action reads: 

S{Q} = II StrQ{i)Q{i + 1) + ^ StrQ(i)A, (4) 

i i 

in terms of which the density-density correlator (||) is expressed as: 
/C(ni,n2;£;) =-lime-pCo / l[dfi[Q{t)] (Q^,^,nK)Q'Mi("2) + Cx<5„,,„,Q12 ^^(nOg^i ^.(ni) exp(-5{g}). (5) 

where the cr-model supermatrix field Q satisfyes a constraint: = 1. With use the Efetov's notations the supermatrix 
Q consists of four supermatrix blocks 

where upper indices 1 and 2 correspond to the "retarded" and "advanced" components of the original RBM model. 
The bottom indices bb in eq.(|^) denote boson-boson components inside the superblocks. For more detailed explanation 
see, e.g. reviews pT|,p|. 

The difference between the two symmetry classes is reduced to symmetry-specific definitions of the model parameters 
p, 7, Co, Cx PI as well as of the supermatrices A and Q. 



^Note, that the factor Cx originates from the number of cross-pairing of components 1^2. Since this contribution gives rise 
at ni = 712 only, it is of main importance in inverse participation ratio (IPR). 
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In the orthogonal ensemble the 8x8 supermatrices Q belong to the coset space UOSP(2,2/2,2)/ 
UOSP(2/2)xUOSP(2/2). Their explicit parametrization, including invariant measure diJ,[Q] can be found, e.g. in 
[p] , p^ . The diagonal matrix A is A = diag(l, 1, 1, 1, —1, —1, —1, —1). The density of states p obeys the semicircle law 



P T 



^{El,, E^f/\ \E\ < E^a. - (2Jo)^/^ (7) 
and the "length scale" parameter 7 is jl^ 

7 = = iB(2Jo - E% (8) 

where the band profile parameters are defined as 

00 oc 

■h= ^= E (9) 

r— — oo r — — 00 

The correlator parameters are Co = 2 and Cx 2, correspondingly.. 

In the unitary ensemble the 4x4 supermatrices Q belong to the coset space U(l,l/2)/U(l/l)xU(l/l), and their 
explicit parametrization can be found, e.g. in p^,p[. In this case the matrix A is A = diag(l, 1, —1, —1). The density 
of states p as well as the parameter 7 are given by |^ : 

P - 7^{Ela. - E^f'\ \E\ < = (4Jo)l/^ (10) 

7 = (7rpJo)'S = ^B{AJo - E^), (11) 

and differ from those in the orthogonal ensemble. It is easy to see, that in the unitary case the "length scale" 7 (at 
-E ~ 0) is twice as large and the correlator parameters are in this case: Co = 1, Cx =1. 

At this level of consideration two symmetry classes differ mainly by the rescaling of the length scale 7 and by 
explicit form of the supermatrices Q. 

At the end of this section let us consider briefly the accuracy of the BRM model mapping onto the supersymmetrical 
cr-niodel. The only approximation done at this step is the integration over so-called "massive" modes performed by 
saddle-point method in Gaussian approximation(for more details see, e.g. |pT|,p|). While the separation between the 
saddle points is of order (see, e.g. |pl],p[) 

A ^ npJo ^ JI'\eI,, - E^)lEl,,,, (12) 



the variances of the massive modes are {6PY « Jo/b [||. Hence, the ratio of the amplitude of the massive modes to 
the sadle points separation is of order of 0(6^/^): 



/ X 1/2 

({SPY) 



E^ f E 



6V2(ij2__ ^2)1/2 \b{E,^a.-\E 



(13) 



which, probably, defines the accuracy of the approxiamtion.^ Putting r ^ 1, one may also conclude, that semicircle 
law for the density of states is valid as long as Ai? = {E^ax — \E\) ^ 1/6. 



*Since the number of saddle is, in fact, even, the cancellation of terms linear in r can occur. Then the accuracy will be better, 
of order of 0{b~^). More careful consideration is needed to arrive at definite conclusion. 
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III. TRANSFER MATRIX SOLUTION: FROM THE DISCRETE MODEL TO THE CONTINUOUS ONE. 



Starting with the pioneer work [|lO| , the transfer matrix technique looks as the most adequate tool to solve the 
problem. Forthcoming works [p^ , p-2[ make some detailes be more transparent. The continuous limit version 
reduces the problem to solution of a partial differential equation. As we check below, this equation is the same for 
both the orthogonal and unitary ensemble, with accuracy 0(7"^). The corresponding boundary conditions are also 
the same when calculating the density-density correlator, but can be symmetry dependent in more general case. To 
demonstrate this let us outline basic steps of the approach. 

Following PI consider the N-site one-dimensional chain of matrices {Q{k)}, with a nearest neighbour interaction 
given by eq.fl). For the sake of definiteness we assume that ni < n2 in eq.dq). Having in mind that the integrations 
over different Q{k), k = 1, ...(ni — 1) in eq.(||) can be performed subsequently one after another with the same kernel 

L{Q, Q') = exp{^StrgQ' - ^eStrQ'A}, (14) 
it is very convenient to consider the quantity Y^^\Q; k) satisfying the recurrence relation 

r(i)(Q; fc) = / dfi{Q')LiQ, Q')Y^^\Q'; k ~ 1). (15) 



with the initial condition y(^)(Q; 0) = 1. In general, the integration over Q is quite nontrivial. A great simplification 
occurs in the limit e — > 0, which is exactly the case we are interested in (see the eq.(|^)). As it was clearly demonstrated 
in [ p"4yi^ ] with Efetov parametrization of matrices Q, after taking such a limit the function Y^^\Q\ k) depends only 
on one of its bosonic "eigenvalues" (say, Ai) of the superblock Q^^, which takes values Ai ^ e^^. Although the kernel 
( p^ has a very extensive grassmannian part, in the limit e — > only the term of the maximal grassmanian order 
survives due to the Parisi-Sourlas-Efetov-Wegner (PSEW) theorem^. This is the term of the eighth or fourth order 
in the orthogonal and unitary case, respectively. Perfoming explicit integration over all variables but Ai and using, 
that Y^^''{Q; k) depends only on z = 2eAi, one can rewrite the recurrence relation ( [Tsl ) in a more simple form: 



y(i)(z;fc)= / dz'L^(4)exp(-z)r(i)(z';fc-l). (16) 







The explicit form of the kernel L-y(p-) can be found, e.g. in |T^] and |lj] for the orthogonal and unitary case, 
respectively. 

This recurrence relation can be used iteratively to construct the function Y^^^z] k) up to the point k = ni. At this 
point the integrand in eq.(^) has an extra factor This factor already depends on the whole set of the ordinary 

parameters of matrix Q and on a half subset of grassmann variables (the two a, a* in the unitary case, and the four 
Q!i, a^, Q!2, a2 in the orthogonal case, in notations of papers p^ , p^ ) ^. Omitting the detailes ( some of them can 
be found in |p2|,p[) we have found that in the limit e — > the only nonvanishing contribution to the density-density 
correlator coming from the factor is the term of the maximal order in grassmann variables. It is of the kind 

c(z, ...)Gmax{a.), with 

„ (' ^ — / oiia\a2a2, orthogonal case, . _^ 

G™a:r(aj = j Q,^*^ unitary case, ^ 

and the ordinary coefficient c{z^ ...), which (after averaging over the angle variables in the orthogonal case) is the same 
for both symmetry classes: c{z, ...) (— zz/(2e)) . 

To proceed further, beginning from the point k — rii we start iterations with a new function Yq (Q; k, ni) satisfying 
initial condition Y^\Q;k — 711,711) ~ zY^^^{z,ni)Gmax{ct). The subscript G denotes here, that this new function 



aquires a nontrivial grassmann part. As a result, one could suspect that the recurrence relation (16) is not valid any 
longer and one should return to the more complex relation (|l^. Fortunately, the closer inspection shows that to the 
leading order in e — > the recurrence relation (15) acts on functions Yg^' (Q; /c, 711) — y^^^ (z; A;, 711) • Gmaxia) in such 



^See, for details [|T|,p|. 

^Analogous factor Qtt.ii) with replacing P ~* a, arises at = n2. 
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a way as if the grassmannian factor remains unchanged, while its ordinary coefhcient Y^^"^ {z; k,ni) is iterated by the 

relation (16) with the same kernel, as before. 

The functions Y^^^{z; k) and Y'^'^\z; fc,ni) allow one to calculate the correlator (^). Let us first consider the "direct 

pairing" contribution to eq.(^) coining from the term ii('^)'3fcb ll(0• 
/Cd(nl, ^2; i;) = - lim -epCo / dp.{Q)Y^^\z = 2e\i-N -n2)Qll^^9,nax{a) {-^)Y^^\z = 2e\i;n2,ni)e-\ (18) 

where the only explicit integration over Q left is that done at the site k = n2- Using the PSEW theorem, performing 
all the integration^Jxcept for that over z — 2eAi, and finally taking the limit e — > we come to the expression 

]Cf{ni,n2]E)^pj ^Y^^\z;N -n2) zY^^\z;n2,ni)e~\ (19) 



Note, that this expression is the same for both symmetry classes, but the functions Y'^^^{z; 712) and Y'^^\z; n2, ni) are 
calculated with kernels L^, which, generally speaking, are dependent on the symmetry class. 

Another contribution to the correlator (||)comes from the cross-pairing term nQll n- It is even more easy to 
calculate. One can verify, that at the point n2 — ni the term Qb6.11Qb6.11 ~ '9b6.ii*3b6.ii ^o the leading order in 
Ai ~ 1/e — > 00. Therefore the full correlator can be written as 

lCE{ni,n2]E) ^ K!'^\ni,n2]E) + C^Smnzl^^E^ i^i, rii, E). (20) 

The final simplification arises if one considers the bandwidth b of the matrices to be large enough, so that the 
parameter 7 ~ &^ 3> 1. Let us note, that this condition was actually used in the previous section in order to map the 
original BRM problem onto the sypersymmetrical ct- model (see discussion in (P). 

First of all, the kernels Lj{z/z') obtained in for the orthogonal case as well as that for the unitary case jlj] 
coincide with each other within the accuracy 0(7^^). Then one can pass to the continuous limit by using that the 



kernel has a sharp peak. It amounts to replace the integral recurrence equation by the differential one |13 



^^^^^(^'-)-W(^)(,,.), R.\y^^-y\, . = 1,2; (21) 



dr ' ' I 

where we introduced a continuous variable x — k/2j, a scaled variable y — 2"fz, and the functions 

W^^\y; x) = y(i)(z = y/27; k = 2jx), (22) 

(y. a;, a;i) = 27 • F (2) ^ y/27; k = 27a;, m = 27x1); (23) 

satisfying the initial conditions 

W^«(y;0) = l, W^^\y;xi,xi)^yW<^''>iy;xi). (24) 



Let us stress, that the differential equation (|2l| ) as well as the initial conditions (g4|) are exactly the same for both 
symmetry classesn. 

The correlator (|l8|) in continuous limit takes a form 



4')(:E2,xi;£;) = f ^W^'\y;xN-X2)yW^'\y;x2,x,)e-y^'\ (25) 



2jJ y2 

where we introduced the diniensionless matrix size as xjv = N/2^. 



''Note, that this integrations give the factor 1 and 2 in the orthogonal and unitary cases, respectively. With the factor Co this 
gives a factor 2 for the both symmetry cases. 

®Note, that in the course of calculations of other quantities, e.g. the probability distribution for the density-density correlator 
or inverse participation ratio, the corresponding initial conditions may be symmetry dependent 0. 
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IV. CALCULATION OF THE CORRELATOR ICsiXi, X2; E) FOR INFINITE BANDED MATRICES. 

The way to obtain a solution of eq.(p2[) was in details considered in 11^,0] ■ For the function W^-^^y,x,xo) one has 



(26) 



where Ka{t) is the MacDonald function. In this work we consider the case of the band matrices of the infinite size: 
iV/27 — > 00. In contrast to |^,|| the correlator arguments xi, X2 are taken at the finite separation {x2 — xi), but both 
far from the ends x = xn = N/2"f. In this limit we have 



M^(i)(y; ^1), W^'^ (y; xn - X2) ^ 2y^'^Ki{2v^/^) 



(27) 



since both xi — > oo, and [xn — X2) 00. Constructing the function W^'^'^ (y^ x, xi) with the initial condition ^4 
W<^^\y,xuxi) ^y2y^/'^Ki{2y^/^) we find 



,1/2 



W^''\y-xi,X2)^^ / d:^2J^2sinh(^i.2)e —^''^-^^'^ K,,,{2y^/^) / dttK,,,{t). 



(28) 



Substituting the functions W^^^ and W^'^^ into the correlator expression ( p5| ) one finally arrives at 

00 

-K.f{x2,x,-E) f dyy-^e-* -22/1/2^1(2^1/2) x 



,1/2 



/ di^2i^2sinh(7rz.2)e r^'^^-^^'i^,,, (2yi/2) / dttK,,,{t) 



(29) 



27r^7 



(ij^2 ^^2 sinh(7ri^2) e 



(2;2~2;i) 



dttKi{t)K,^Jt) 



dv V siim(Tny) e " 



16-27 



1 + cosh TTf 



(30) 



In eq.(|29D we omitted the factor exp(— y/27) since the integral over y converges at y 1 but 27 ~ 6^ ^ i •yyg see 
also, that in the infinite matrix size limit the correlator depends on the difference x — X2 ~ xi only which means the 
restoration of the translational symmetry. 

It is interesting, that the BRM model prediction for the functional form of the correlator )Ce{x2 — xi; E) coinsides 
exactly with that obtained many years ago by Gogolin ||] for the strictly ID electron system in disordered potential. 
However, the spatial scales in these systems are quite different. Namely, in the strictly ID electron system it is equal 
to the mean free path defined by the Born amplitude of electron-impurity scattering, while in the quasi-lT) BRM 
model it is governed by the localisation length Ic — 2j (x , growing quadratically with matrix bandwidth. 



V. REM MODEL PREDICTIONS FOR TIME- AVERAGE CORRELATOR ICt{Xi,X2). 

As we have already mentioned in the Introduction, the time-averaged correlator ICt{xi,X2) can be obtained from 
the "fixed energy" correlator ICe{xi,X2; E) by integrating the latter over the energy: 



00 

K.t{xi,X2)= J dEK.Eixi,X2;E). 



(31) 



Substituting eq.(30) into eq.(|3l|) we get 
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lC,{x,-x,)= j d..smh(,r.).exp| ^ ,,)| j , (32) 

-So 

with energy dependent localization length lc{E) = 2'-f{E). The scaled variables xi,2 are related here to the discrete 
ones ni^2 via the energy independent scale factor lc{0)- xi^2 — '^i,2/'c(0). 

For large N,b:^ 1 the density of states p{E) in RBM model obeys the semicircle law: 

piE) = piO) . (1 - I) - 4 (1 - I) / dE piE) . 1; (33) 

and the localization length is, correspondingly, 

UE) = UO) • (1 - ^) (34) 

with Eq = 2Jo, ^0(0) = JqB/2 « or Eq = 4Jo, /o(0) = BJq w 26^/3 for the orthogonal or unitary cases, 

respectively (see eqs.(^-(^). 

Let us make a remark about the limits of the integration. According to the note at the end of Sec.||, both the 
semicircle law for the density of states and the solution itself are valid as long as {Eq — \E\) 3> SE ~ Eo/b. Thus, 
extending the integration in the Eq.(^2|) over the whole region [—Eq, Eq] we get an uncertainty coming from the regions 
Eq — \E\ < SE, where our integrand differs substantially from the exact one. To estimate an order of magnitude 
of this uncertainty, let us consider the contribution coming from this regions with our integrand. In the "worst" 
case of the coinciding correlator arguments xi = xi the integrand has integrable square root singularities at the ends 
E ~ zLEmax, and their contribution has a relative order of mag nitude {SE/EqY/^ - in the case of nonzero 

X = X2 — xi the contribution of of the regions Eq — \E\ < 6E is exponentially small at x ^ b^^. In particular, we 
conclude that in the recent numerical BRM simulations done for the matrix bandwidth b — 4—12 substantial 
deviations order of b~-^/^ can occur for the correlator at a; ^ b^^. 

The integration in eq. (|3^ ) can hardly be performed analytically. Let us therefore consider two simple limiting cases 
first. It can be most easily done for the correlator value at a; = (12 — 0:1) = 0: 

/Ct(0) = 2/CB(0;£;)b=o=^^(l + Cx), (35) 

which is a "time-averaged" version of the so-called inverse participation ratio (IPR) £,t = I'/'al^- Note, that in the 
BRM model "time-averaged" IPR is twice as large as compared to its "fixed energy" counterpart usually considered in 
the solid physics and defined as = X]q Ifal'^ S{Ea — E) for quantum states near E ~ 0. It is a natural consequence 
of the averaging over localized quantum states with lc{E) < Iq{0). 

Large-x asymptotics {x ^ ^c(O)) can be calculated by the steepest-descent method. Performing the calculation, we 
find 



K-tix) 



4 //c(0)V/' 



K.e{x;E) |_e=o ■ (36) 



Comparing with the expression eq.(^5|)we see that at large x the extra factor reduces the value of the correlator. 
Therefore, converting K-e into /Cf changes the correlator form, most considerably in the region ^ ^c(O). 

In the recent numerical work the data for the correlator K,t{x) have been obtained in the course of a direct 
numerical simulation of the BRM orthogonal ensemble. With use of the Gogolin formula |^ the best fit value for the 
localisation length was found Ic ~ 0.296^. In Fig.l we plot the correlator lCt{x) predicted by eq.(|3^) which has no 
free parameters and corresponds to lc{E = {)) = b"^ /"i ^. The dashed line presents the parametrization of the data by 
the Gogolin formula with Ic — 0.296^. To illustrate an accuracy of the agreement, we plot also the Gogolin formula 
with Ic — b^ /'i. We see, that the agreement looks reasonable. For more reliable analysis one needs to compare the 
correlator JCt{x) directly with the data for BRM ensembles with the matrix band with b':^ 1. 



^To facilitate the comparison with data we plot the curves in terms of scaled variable x — nj}? , as defined in Ref. 



7 



However, the most crucial test of our theory would be the comparison of the correlator behavior at a; < 1 with 
numerical data (see, Fig. 2). We expect it to be essentially more peaked at a; = than the "fixed energy" correlator 
K-e{x2 — xi]E) given by the Gogolin formula with Ic = \? jZ. Apart from this fact, it would also be interesting to 
compare the time-averaged IPR with data, which is predicted to be twice as large in comparison with the corresponding 
"fixed energy" quantity. 



VI. CONCLUSIONS. 



In this work we have studied the density-density correlators in the BRM ensembles for the orthogonal and unitary 
symmetry cases. In contrast to the works [p|j5|,|9|,|l3| we consider the case of infinite matrix size N ^ oo, with correlator 
arguments ni < n2 taken well inside the argument region: [1, ni, n2, iV], far from the ends: ni S> b^, 
TV — 712 ^ 6^ {b is matrix bandwidth). 

We have found, that the "fixed energy" density-density correlator (|^) turns out to be the same in both the orthogonal 
and unitary BRM ensembles, within the accuracy 0{b~^). Moreover, its functional form coincides exactly with the 
old Gogolin formula originally derived for ID electron system in disordered potential. 

It deserves to be mentioned that the analytical solution of the BRM model given in consists of several approx- 
imation steps. The first one is the mapping of the model onto the supersymmetrical cr-modcl. It has, probably, the 
worst accuracy of order 0{b~^^^). At the same time the accuracy of the cr-model solution, including the continuous 
approximation, is the order 0(7"^) ^ 0{b^^). Thus, one may expect deviations from the exact solution to be of 
order ofO{b-^/^). 

Beside the "fixed energy" density-density correlator JCe{x2 — Xi;E), we have also calculated the time-averaged 
correlator /Ct(x2 — xi). Two functions are found to be different, especially in the region of {x2 — xi) < 1. We predict 
also, that the time-averaged inverse participation ratio for BRM model is twice as large as the "fixed energy" one. 
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FIG. 1. The time-averaged density-density correlator JCt{x) versus x — n/lP . The solid line represents the prediction by 
eq.(^). The dashed line shows the parametrization of the numerical data[3], the Gogolin formula with Ic = 0.29, and dotted 
line corresponds the same Gogolin formula with Zc = 0.33. 
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FIG. 2. The time-averaged density-density correlator K,t{x) versus x — nj}? (central region). The solid line represents 
prediction by eq.(p2|). The dotted line shows the Gogolin formula with 1^ = 0.33. 
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